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We consider the thermal correlation functions of vector and axial-vector currents and evaluate 
corrections to the vector and axial-vector meson pole terms to one loop in chiral perturbation 
theory. As expected, the pole positions do not shift to leading order in temperature. But the 
residues decrease with temperature. 
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O I. INTRODUCTION 

o . 

(N ■ 

An important topic in strong interaction at non-zero temperature is the calculation of shifts in the masses and the 
couplings which hadrons suffer while moving through a heat bath. Apart from their immediate relevance in under- 
standing the experimental data on heavy ion collisions, these shifts as functions of temperature provide an effective 
description of the thermal properties of strong interaction. Mention must be made here also of the condensates, which 
are related more directly to the condensed matter aspects of strong interaction, particularly to the possibility of phase 
transition with the rise of temperature. 

At low energies processes involving strong interaction as described by QCD are best analysed in chiral perturbation 
theory It has long been established phenomenologically that the symmetry group SU(3)r x SU(3)l of the QCD 
| Lagrangian of three flavours of massless quarks is broken spontaneously to its diagonal subgroup SU (3)y, giving rise 
. to the octet of Goldstone bosons. Chiral perturbation theory realises this symmetry nonlinearly on the physical fields. 
£f) ' This framework restricts considerably the interactions among hadrons involving the Goldstone bosons. It has led to 
, a whole set of accurate analysis of low energy hadronic processes . 

Chiral perturbation theory finds a natural application to problems of strong interaction at non-zero temperature 
T. At low temperature pions dominate the heat bath, the other hadrons being exponentially suppressed. One is thus 
led to consider the reduced chiral symmetry group G = SU(2)r x SU(2)l, which breaks spontaneously to the isospin 
Oh| subgroup H = SU(2)y. This symmetry implies that interactions involving low energy pions are weak. In calculating 
quantities at low temperature these interactions can therefore be treated as perturbation. A wide variety of problems 
has been treated in this way pl[|. In particular, it has been applied to the correlation function of nucleon currents 
to find the position of the nucleon pole and its residue at finite temperature 0]. In this work we apply a similar 
procedure to the isotriplets of vector and axial-vector mesons. 

Here we consider the correlation functions of vector and axial- vector currents defined in QCD. To calculate the 
pole parameters to 0(T 2 ), we draw all the one loop Feynman diagrams with vector and axial- vector mesons and 
JH I pions, which contribute to these correlation functions. With the form of interaction vertices prescribed by chiral 
perturbation theory, it is simple to evaluate the finite temperature parts of these diagrams near the poles in the chiral 
limit. Although we begin with a large number of diagrams, only a few actually contribute to the meson poles to order 
T 2 . We also argue that other particles not included in our scheme cannot change the results to this order. 

In sec. II we summarise the chiral transformation rules of the relevant field variables to arrive at the chirally 
invariant Lagrangian of the spin one mesons and pions. In sec. Ill we obtain explicitly the interaction vertices needed 
to evaluate the Feynman diagrams. In sec. IV we find the shift in the masses and couplings of these mesons to 0(T 2 ). 
Finally in sec. V we summarise our results and comment on other attempts to determine these shifts. The Appendix 
discusses the 2x2 matrix and the Lorentz tensor structures of the correlation functions in the real time thermal field 
theory. 



II. CHIRAL PERTURBATION THEORY 



It is convenient to review in this section the definition of the appropriate field variables and their transformation 
properties, leading to the effective Lagrangian of chiral perturbation theory The Goldstone fields -K a {x) reside 

in the coset space of G with respect to H . A standard parametrization of this space allows one to collect these fields 
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in the form of an unitary matrix, 

3 

u{x) = e ™M/ 2F " , ttO) =^K a (x)T a , (2.1) 

a=l 

where the constant F„ can be identified with the pion decay constant, F v = 93 MeV and r a are the Pauli matrices. 
The matrix u transforms under G according to the rule, 

u -> g R uh^ = hug\ , (2.2) 

where g Ri L £ SU(2) Ri l and h(n) £ SU(2)v- The matrix U — u 2 then transforms as 

U^g R Ugl (2.3) 

Any multiplet i/j(x) of non-Goldstone fields transforms as 

ip -> hip, (2.4) 

with h in the appropriate representation. In particular, if they also belong to the adjoint (triplet) representation, such 
as the vector and the axial- vector mesons we are considering here, we may use the same h as in (2.2): Denoting the 
triplet fields by i? M (x), 

i 3 

R M (x) = -j=52l%(x)T», (2.5) 

V a— 1 

it transforms as 

Rp, — * h h) . (2.6) 
The singlet field S^(x), of course, remains unchanged, 

As already stated we are concerned here with the evaluation of the two-point functions of vector and axial-vector 
currents of QCD, 

V?(x) = q(x)^q(x), Aftx) = q(xh^q(x), (2.7) 

in the effective field theory. This is most conveniently carried out in the external field method, where one introduces 
external fields v^(x) and a°(i) coupled to the currents V£(x) and A®(x) (l). The original QCD Lagrangian Cqq D of 
massless quarks is then augmented to 

C$l D +Wx)V£(x)+a*(x)AZ(x) 
= iqRj^id^ - i{v^ + a^)}q R + iqL-f^id^ - i(v^ - a AI )}q L H (2.8) 

where q R {qL) is the right (left) handed component of q, q R} L — s(l ± 7s)<7- The ellipsis denote the (flavour neutral) 
gauge field term, and v^(x) and a^(x) are matrices in flavour space, 

^W=<Wy . a M (a:)=oJ(i)y. (2.9) 

The global invariance of £qp£> is now raised to local invariance under the transformations 

Qr = 9R9R . 1l = 9LqL, 
v 'n ± °V = 9r,l( v h ± a »)9 R , L + i9R,Ldf,g R L , (2.10) 
where the group elements g Rl h are now x-dependent. 
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The field strengths corresponding to the external potentials are given as usual by 

F^ L = a"(«" ± a") - ± a") - i[/ ± / ± a"]. (2.11) 
The covariant derivatives of U and i? M are given by 

D^U = d^U - i( Vfl + ap)U + iU(v^ - a„), (2.12) 

and 

V l Jt v = d vl B v + \C^R u ] (2.13) 

with the connection 

r '' = \ - i<yVfJ * + a ^ u + _ ^ Vfi ~ a ^)] wt ) ( 2 - 14 ) 

Thus the building elements of the effective Lagrangian are U, D^U, F^ L , and V^i?^. 

Note that while some of the above variables (U, D^U, F^ L ) transform under the full group G, others and 
V M i?„) transform under the unbroken subgroup H. One may take advantage of the mixed transformation rule of u 
to redefine the former variables so as to transform under H only. Thus one introduces the field variables |7j], 

Uf, = in* D flu 1 = ul (2.15) 

and 

/£" = ±u ] F» v u + uF£ v v) (2.16) 

We now use these variables to write down the leading terms of the different pieces of the Lagrangian density, which 
are invariant under H and hence also under G. To calculate the correlation functions of isotriplets of vector and 
axial- vector currents, we need the Lagrangian for the interacting system of both the isotriplets and isosinglets of 
vector and axial-vector mesons and pions in presence of the external fields. (We shall see later that other particles 
cannot contribute to the pole terms to order T 2 .) For the pions we have the familiar term, 

£M = <«„«">• ( 2 - 17 ) 

Here and below the symbol (A) stands for the trace of the 2x2 matrix A. For the spin-1 meson isotriplet and isosinglet 
fields, let 

Rpv = V p,R v — Vy-R^ , Spv — dfl v — d v S,j,. (2-18) 
Then the kinetic terms take the form 

C km {R, S)=Y, (-\(R^Rn + lm%(R p Rn) + £ (- \ + \m 2 s Sfl^ , (2.19) 

where the two sums run over the isotriplets and the isosinglets respectively. 

The leading coupling terms linear in the octets of the vector and the axial-vector fields have been obtained in Ref. 
jjlD for the symmetry group SU(3) x SU(3). As mentioned already, the symmetry reduces to SU(2) x SU{2) at finite 
temperature. In effecting this reduction, we keep the masses of the physical particles in each of the octets to be the 
same but consider the interaction separately for the isotriplets [p(770), ai(1230)] and isosinglets [a>(782), /i(1282)]. 
Restricting to terms relevant to one loop calculations, we get for the vector meson couplings, 

1 

2V2tov ' 

while for the axial vector meson such couplings are 

£int(A) = _A_(F 01 <a 1/t ,,/r> +iH ai (a llM [u u J^]) + H h e^ Xa f?{u v f>l°)). (2.21) 

Finally we write down the quadratic couplings of the triplets with the singlets and between themselves, 

C mt (V,A) = ^e^gru^p^u^+gift^u'*)) + l -g 3 {p^[a lix ,u v ]). ( 2 - 22 ) 

(The term (a^ [pn,u v ]) is equivalent to the third term above to leading order in pion momentum.) Although we 
have here a number of different coupling constants, only some of them will actually enter our results for the shifts in 
particle parameters to 0(T 2 ). 



£int(V) = — {F p { Pllv f» v ) + iG p (p flv [u fl ,u u }) +iH p { Pli [u v J» v ]) +H u e^ x ^(u v fl a )), (2.20) 
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III. PERTURB ATIVE VERTICES 



The vertices needed in evaluating the one loop Feynman diagrams for the correlation functions can be recognised 
best from the diagrams themselves. Let us consider the correlation function of two vector currents at non-zero 
temperature T = 1/(3, 



I 



d i xe lqx TrgTV^{x)V^(0) 



(3.1) 



where g — e^^ H /Tre f3H is the thermal density matrix of QCD. The diagrams contributing to it are of three types, 
namely those of self-energy, vertex modification and intermediate states shown respectively in Figs. 1, 2 and 3. Wc 
see that we only need vertices up to four fields, counting both quantum (particle) and classical (external) fields. At 
each vertex the pion field can appear at most quadratically. When it does appear so, the two pion fields must be 
contracted at the same vertex. The resulting pion loop is of 0(T 2 ), if no derivative is present on the pion fields; but 
the loop is of higher order if it does, allowing us to ignore such vertices. Also at each vertex the external fields can 
occur at most linearly (with the exception of the vertex in Fig. 3(b)). At such a vertex with an external field, the 
vector or the axial vector meson field must also occur linearly. Keeping these requirements in mind, we now obtain 
the polynomial version of the chiral Lagrangian written in the previous section. 



<*> Wl/V 



(0, a 1( 7C 



«> MXPJV 




WW 



FIG. 1. p meson pole and self-energy diagrams 
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FIG. 2. Vertex correction diagrams 
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(b) 




FIG. 3. Intermediate state diagrams 



The pion Lagrangian (2.17) can be expanded as 



+ ^F 2 a^ ■ a* + i {tt • 7r(w^ ■ - o M ■ a**) - u M • nv^ • tt + a M • ttc^ • tt} . 



(3.2) 



Here and below we denote the isovector triplet of fields by bold face letters. Note that we have added the pion mass 
term, even though we shall evaluate the final results in the chiral limit. The first term gives the free, thermal pion 
propagator matrix. As discussed in the Appendix, all we need here is its 11-component, 



d 4 xe lk - x TrgT7T a {x)TT b (0)\ 11 



gab 



k 2 — m 2 + ie 



- + 2irn(k a )5(k 2 - ml) 



(3.3) 



where n is the bosonic distribution function n(fco) = (e' 3 ' feo ' — 1) 1 . 
The kinetic part of the Lagrangian for the spin one mesons reduces to 



C kin {R,S) = -\ J^l^R, ■ (d^R" - d v B?) - m 2 R R„ ■ R»} - \ J^S^S" - PS") - m 2 s S^}. (3.4) 



p,ai 



The density of spin one mesons are suppressed exponentially in the heat bath. So we need only the vacuum part of 
the respective 11-components of the propagator matrices, which for J?" for example, is 



/ 



dtxjP* (0|Ti$(a;)i£(0)|0) = S ab l{ 9^+P^Pu/m 2 R ) 



(3.5) 



As an illustration we now work out the vertices present in the term (p^f^)- Since is proportional to the 
external fields, we need only the terms in p^ v without them. Thus upto four field vertices, we have 



(3.6) 



As already explained, we have to contract the two pion fields in the four field vertex above. Its thermal piece is 



TrQTn a (x)* b (x)\u - 5 ab J ^L n (k )6(k 2 ml) = 5 ab ^. 



(3.7) 



in the chiral limit. (We see here that any derivative on the pion fields would result in contributions of higher order in 
T). Thus the third term in Eq.(3.6) reduces to the first with a T-dependent coefficient. 
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as 



In this way we derive the pieces of the polynomial Lagrangian from those of the chiral Lagrangian of Eqs. (2.20-22) 



Ttly 



2a " ^ Pv ■ d»n x ff'ir - ^^-e^xo^d"* ■ 0V, (3.8) 



C{A) = ( ( 1 - 0"a" • (c^ - d„a llt ) + ^-0"t>" ■ (d^a lv - d v a^) x tt) 

+ ^S L ^Xaf^d^-d x a\ (3.9) 



and [ud p v = (d^u)v — ud p v 



C{V,A) = e^xa [ ^d t p x ■ d a n+^f^a x ■d°Tv)+ §-d»p" ■ (a lfM x 9„tt - a Xv x m tt). (3.10) 

In writing the above expressions, when two or more vertices present themselves with the same field content, differing 
only in the number of derivatives on the pion field, we retain only the one with fewer or no derivative on it. 

The coupling constants can generally be determined from the observed decay rates of the particles. Thus the decay 
rate T(p° -> e + e~) = 6.9 keV gives F p = 154 MeV. Similarly the decay rate T{a x -> 7rV) = 640 keV gives F ai = 135 
MeV. The latter constant can also be determined from one of the Weinberg sum rules [[|, agreeing closely with this 
value. The decay width T(p — > 2ir) = 153 MeV gives the coupling G p — 69 MeV. Here one expects large chiral 
corrections as the pions are far from being soft. Thus the chiral loops reduce this value to G p = 53 MeV 0. Using 
this value and the decay rate r(w — * 3ir) = 7.6 MeV, we get the dimensionless coupling constant g\ — 0.87 p0[ . 
There does not appear any data relating to the decay of f± to determine (72 ■ The remaining coupling constants in 
Eqs. (3.8-10) will not appear in our results to 0(T 2 ). 



IV. MASS AND COUPLING SHIFTS 



Given the vertices, it is simple to calculate the diagrams of Figs. 1-3. The Lorentz tensor and the thermal matrix 
structures of the two point functions are discussed in the Appendix, according to which the free pole term of Fig. 1(a) 
is given in the variable E = qo for q = essentially by 



-£ 4 



{F p /m v f 
E 2 — my + it 



We now determine how this pole position and the residue are modified by interactions at finite temperature. The 
self-energy diagrams of Fig. 1 modify it to 

_ e4 {F p /myf 



E 2 - m\ - U t 

where we use Eqs. (A. 12-13) to construct II t from the finite temperature part of the polarisation tensor IIq,^, 

the sum running over the diagrams of Fig. 1(b). (We have left out Fig. 1(c), which is clearly of order T 4 .) To write 
out these contributions we define the gauge invariant tensors 

A a p(q) = -g a p + qaqp/q 2 , 

B a0 (q, k) = q 2 k a kp - q ■ k(q a kp + k a qp) + (q ■ k) 2 g a/3 , 

C a p{q, k) = q 4 k a k/3 - q 2 (q ■ k)(q a k/3 + k a q p ) + (q ■ k) 2 q a qp. (4.2) 
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The three contributions are now given by 



n H (29l\ f d*k n(k )S(k 2 -ml) 2 2 



tt( o nf^Y f n(k ) S(k 2 - mp 2 

a/3 = ( Tt - ) / 7^3 -fZ UY2 Z^T~ ~ C a/3 / 'm A ) , (4.4) 



and 



U a/3 = 2 I ^-^j ) / — M2 _ to2 C Q/3 . (4.5) 



The respective contributions to Il t are {u>k — y k 2 + m 2 ), 

d 3 kn(oJk) k 2 



^\E) = ^E 2 (E 2 - m 2 v + m l)J^ 



f2uj k {E 2 - m\ + ml ) 2 - 4E 2 u 2 ' 



(4.6) 



U^ { E) = U^) - m \ + ml) [ ^ ■ , Jf(* + ^ + , , (4.7) 



3 V-*W V A J (2tt) 3 2w fc (£ 2 - m 2 + m 2 ) 2 - 4E 2 



k 



and 

n * (£) -H^fJ Nw^'rai' (48) 

To evaluate the integrals for E 12 near m 2 / in the chiral limit, we note the difference in the small k behaviour of the 
denominators of the integrands: while it behaves like ~ k 2 in Eq. (4.6), they are constants in Eq. (4.7-8). Accordingly 
we get 0, 

H^{E)=-^(E 2 -m 2 v ), n^~nW~0(T 4 ). (4.9) 
Including also the constant vertex corrections from Figs. 2(a,b), we get for the complete pole term 

-E-^l ,4.10, 



with 



3 J \2F 2 \ 



^ = F P \l-(l + ^)^\. (4.11) 



It is simple to see that the remaining diagrams cannot alter this pole term to 0(T 2 ). Each of the remaining vertex 
corrections of Fig. 2(c,d) is essentially of the form of the corresponding self-energy integral (lit) multiplied by the 
p-meson pole. Then the behaviour (4.9) of the self energies excludes any pole contribution to 0(T 2 ) from these 

corrections. Of the contributions from the diagrams of Fig. 3 with intermediate states, the one with ttuj is like Il\ , 

while those with irai and tttt have two powers of k fewer in the integrands compared to the self energies n| Ql) and 

H t respectively. Clearly they cannot contribute to the p meson pole. 

It is now also simple to establish that any particle with a mass other than my and not included in our calculation, 
cannot alter our results. The argument rests simply on the structure of the interaction vertices without external fields. 
Such a vertex must have a pion field with a derivative. Thus the loop integrals involving such a particle must be of 
order higher than T 2 . 
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It is interesting to work out the amplitude of Fig. 3(a) corresponding to the intermediate state irai. Its finite 
temperature part is given by 

F ai V f ^)^ 2 -^ { f {(f „ 2q . k)A ^„ B ^ (4 . 12) 



m A F n J J (2tt) 3 (q-k) 2 -n 
It may be readily evaluated in the chiral limit to give essentially [O 



T'l^iF^f (4 , 3) 



6F 2 [ E 2 

where the expression in second bracket is the contribution of the axial vector meson a\ to the vacuum correlation 
function of axial vector currents. Observe the different origin of the 0(T 2 ) correction to the p and ai poles. While 
for p, it is given by the self-energy due to ttlu and constant vertex diagram, it is given by the intermediate state wai 
for a\. This presence of axial vector meson pole in the vector current correlation function at finite temperature was 
detected earlier using PCAC and current algebra Q. 

One can calculate the axial-vector meson pole term in the axial-vector current correlation function in an entirely 
analogous manner. Again one finds that the pole position does not shift from E 2 = m\, and the residue is given 
by Eq. (4.11) with g\ replaced by gi. 



V. DISCUSSION 



Here we have determined the temperature dependence of the parameters of p and a\ meson poles appearing in the 
vector and the axial-vector correlation functions. It is based on the flavour symmetry of QCD and its spontaneous 
breakdown, as embodied in chiral perturbation theory. It allows a low temperature expansion of these quantities, 
of which we just calculate the leading term to 0(T 2 ). To this order the mass shift is zero, due to the presence of 
derivatives on pion fields in interaction vertices as required by the chiral symmetry. The residues depend on the 
coupling constants at the vertices pun and a±fiw and decrease with temperature. These results parallel those for the 
nucleon calculated earlier, the above vertices playing roles analogous to that of the vertex irNN Q . 

We may try to estimate the range of validity of our results by comparing these with the temperature expansions of 
other quantities. For the quark condensate the expansion has been worked out to order T 6 p|. It shows that the T 2 
term continues to dominate the series up to T = 150 MeV. It may not be true in all cases, however. Take the example 
of the nucleon pole residue in the nucleon correlation function. In the pion-nucleon system there is the low energy 
resonance A(1237), which limits considerably the range of validity of the first few terms in the chiral perturbation 
series for irN scattering amplitudes. Already at T — 100 MeV, the average pion energies in the heat bath are outside 
this range. This situation limits the validity of the T 2 term in the nucleon pole residue to temperatures below 100 
MeV. 

As the Tip and 7rai systems do not have any such low energy excitations, we may expect our calculation to order 
T 2 to represent the expansion well up to temperatures of about 150 MeV, as in the case of the quark condensate. 
There is, however, one source of enhanced chiral corrections to our results, which may limit this range of validity. 
As mentioned earlier, at the vertices where the p or the external fields couple to two pions (see Figs. 1(b), 2(c,d) 
and 3(a)), the pion momenta are not at all soft. Thus, although these diagrams do not contribute to order T 2 , the 
coefficients of higher order terms may be larger. 

We finally comment briefly on the different approaches in the literature to determine the temperature dependence 
of these parameters. First, there are other Lagrangian approaches. One has the massive Yang-Mills |l3| and the 
hidden gauge |Q formulations. Since these Lagrangians are built on the SU(2) x SU(2) flavor symmetry, they also 
produce no shift in the meson masses to 0{T 2 ) Jl5| , p6[ . However, the coefficients of T 4 quoted in these works must 
be incomplete, as they do not include two loop diagrams. One has also the more phenomenological Lagrangian of 
Quantum Hadrodynamics [jl7|, which produces negligible shift in the meson masses to 0(T 2 ) 

The second approach tries to calculate the correlation functions using soft pion techniques, but fails to reproduce the 
terms proportional to g\ and g\ in F p and F ai respectively [Q. Since chiral perturbation theory actually incorporates 
such techniques in a systematic and refined way, it is interesting to trace this disagreement. The Tr(ace) in Eq. (3.1) 
can be restricted to the vacuum and the one pion states, giving the full contribution to 0(T 2 ) in the chiral limit. In 
the pion matrix element one may use the hypotheses of partially conserved axial-vector current and the algebra of 
currents to arrive at 

TrgTV2(x)V v b (0) = (l-^\ (0|T V? (x)V^ (0)|0> p + ^(0^^(^(0)10), 
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and a similar one for the axial vector correlation function. The subscript [3 in the first matrix element on the right 
is a reminder that it is really not a vacuum matrix element, since any internal (pion) line which would occur in its 
perturbative evaluation must be taken as thermal. Being multiplied with T 2 , the second matrix element can be taken 
as a true vacuum expectation value for results to 0(T 2 ) and we have already verified the correctness of this term in 
our calculations above. 

The authors of Ref. |]l2| assume the first matrix element also to be a true vacuum matrix element having no 
temperature dependence. In terms of Fcynman diagrams it amounts to ignoring the self energy diagrams of Fig. 1, 
which can in general contribute to both the pole position and the residue. In fact, the missing g\ term in F p arises 
exactly from the self energy diagram of Fig. 1(b). A similar situation arises also in the case of the correlation function 
of nucleon currents. Applying the same soft pion techniques and assuming the resulting matrix elements to be true 
vacuum expectation values, one misses the term proportional to the axial coupling constant of the nucleon in the 
residue of the nucleon pole [^9[ . We conclude that such soft pion techniques as applied to the correlation functions 
cannot reproduce the full results of chiral perturbation theon 



The third approach is based on thermal QCD sum rules pO 21. It is convenient to discuss these sum rules after 
subtracting out the corresponding vacuum sum rules, leaving only T dependent contributions. The spectral side of 
such a sum rule is given by three types of one loop diagrams, like those of Fig.l, 2 and 3. In the literature, however, 
it is assumed to be saturated by a pole term with T dependent parameters and the two particle intermediate states 
containing (at least) one pion. The inadequacy of this saturation scheme is revealed by the fact that the self-energy 
and vertex correction diagrams are generally not exhausted by contributions which make the pole term T dependent. 
It is the neglect of the remaining contributions which make the existing results unreliable. (For incorporation of such 
remaining contributions in the nucleon sum rule, see Refs. Jl9| , p2] |.) 

The fourth and the last one is a phenomenological approach based on the (first order) virial expansion of the self 
energy of a particle |],§3|]. Quite generally it expresses the complex shift in the pole position of a particle as an 
integral over the product of the forward scattering amplitude of a pion on the particle and the pion distribution 
function. Since the pion gas is rather dilute even at a temperature ~ 100 MeV, this formula should hold upto higher 
temperatures than the results of chiral perturbation expansion to first order only. But the problem here is to construct 
the scattering amplitude, which is reliable p4| ]. 

APPENDIX 

Here we discuss the kinematic structure of the vector and axial-vector two point functions. In the real time 
formulation of thermal field theory we are using here, one has with each physical vertex an associated (ghost) vertex 
to take into account |2q] . This causes all two point functions to take the form of 2x2 matrices. Thus even if we begin 
with the 11-component, the perturbative expansion will mix it up with the other components. But as we show below, 
this complication does not arise for the problem at hand and we actually need work with the 11-component only. 

Consider first the complete vector (or axial-vector) meson propagator, which appears in the self-energy diagrams 
of Fig.l. Denoting the thermal 2x2 matrices in this Appendix by bold face letters, it satisfies the Dyson equation, 

IV = D$ +Dfl(-iU x nD <ru (A.l) 
The free thermal propagator has the factorised form 



D$(Qo,q) = U( qo ) ( D ® J oU ) U(q ) (A.2) 

where 



D (0): 



Ul ^ = i^ Jr^) • " 1 1 (A ' ; - !1 



/ n 



and DjxJ (q) is the vacuum propagator 



pm) = i{ -~ 9 T + Y v,m2R) (a.4) 

M q z — m R + ie 

From its spectral representation, the complete propagator iD^^q) can also be shown to admit a similar factorisation, 

D^(q ,q) = U(q ) ( D g" £ )u(q ) (A.5) 







It then follows from Eq. (Al) that U(— iTl^U must have the diagonal structure, 



\1V 

The matrix equation (A.l) now collapses to the ordinary equation 



U(-iU au )U = ( ^ ^ ) (A.6) 



D^^DW+D^-m^D^ (A.7) 

and its complex conjugate. The only remnants of the matrix structure are the relations between and the 
components of Tl uv . With the 11-component of the latter, these are 

Re(ILu) uv = ReU^ , 7m(n n )^ = (1 + 2n)ImU^. (A.8) 

In calculations at finite temperature, one usually prefers the matter rest frame, thereby losing explicit Lorentz 
covariance. The latter may be restored by bringing in the four velocity in the medium Q Then the time and space 
components of a four-vector q M are converted to Lorentz scalars, ui = u ■ q and q — \/oj 2 — q 2 . In this framework a 
gauge covariant decomposition of the polarisation tensor is 

IL li „{q) = n t P„ u +IL l Q IM/ , (A.9) 

where we choose the kinematic covariants as 

Pan = ~9imu H n ^Uullv , Quu — ~Z^UhU v , (A. 10) 

q q q 

with — u^ — bjq^/q 2 . The invariant amplitudes are functions of two variables, say g 2 and to. With the decomposition 
(A.9) the Dyson equation (A.7) can be solved to get 

n - I p i \ Qjuj_ , _J_ Qjj>Qv_ 

q A — m R — Lit + ie Q ttt-r — 9 Hi + «e ? ? 

To find the invariant amplitudes IT t) / we form the scalars 

n 1 = n£ , n 2 = uW!i uv , (A.12) 

when Eq. (A.9) gives 



n, = 4n 2 , u t = -hnx + ^n 2 ) . (A.13) 

<7 Z 2 q 



In the limit q — > 0,the kinematic structures (A. 10 ) relate the two invariant amplitudes, 

n t (g ,<f= 0) = q 2 TL l (q 0l q = 0) . (A. 14) 

In this limit the second equation of (A.13) simplifies to 

n t = -in. (A.15) 

Having discussed the 2x2 matrix structure of the complete propagator, we turn to the same for the two point 
functions, which for definiteness we take the one for vector currents. From the self-energy diagrams of Fig. 1 it gets 
the contribution 

S ab (J^\ (q 2 9^ - q^q X ) D Xa (q 2 g° v - q° q v ) 

= S ab U(q ) ( G ^ _jL ) U(q Q ) (A.16) 



No confusion should arise from the earlier use of as a field variable in sec. II 
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where 



my 

on using (A. 11). Note that at q = the two pole positions above coincide due to the kinematic constraint (A. 14). 
The vertex correction diagrams of Figs. 2(c) and (d) are of the form 

nMB S = ( ^f 2 n J D w* ) (A.18) 

while the diagrams Fig. 3(a) with intermediate states are of the form given by (A. 5). 

We see that the results of evaluation of different diagrams are typically of the form of a diagonal matrix, sandwiched 
between U or U~ . Near the meson pole the U matrix reduces to the unit matrix with exponential correction 
~ e~ m p/ T '. Such corrections are too small to be retained when we are calculating corrections of 0(T 2 ) only. The 
remaining diagonal matrix may be represented by its 11-component, which again admits a decomposition into invariant 
amplitudes like that of in Eq. (A. 9) with the longitudinal and the transverse ones being again related for q = 
as in Eq. (A. 14) for n /it . 
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